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ABSTRACT 

We present a spectral investigation of X-ray binaries in NGC 5128 (Cen A), using six 100 ks Chandra 
observations taken over two months in 2007. We divide our sample into thermally and non-thermally dominated 
states based on the behavior of the fitted absorption column Nh, and present the spectral parameters of sources 
with Li ;J 2 X lO-'^ erg s"'. The majority of sources are consistent with being neutron star low mass X-ray 
binaries (NS LMXBs) and we identify three transient black hole (BH) LMXB candidates coincident with the 
dust lane, which is the remnant of a small late-type galaxy. Our results also provide tentative support for the 
apparent 'gap' in the mass distribution of compact objects between ~ 2-5 Mq. 

We propose that BH LMXBs are preferentially found in the dust lane, and suggest this is because of the 
younger stellar population. The majority 70-80%) of potential Roche-lobe filling donors in the Cen A 
halo are > 12 Gyr old, while BH LMXBs require donors > IMq to produce the observed peak luminosities. 
This requirement for more massive donors may also explain recent results that claim a steepening of the X-ray 
luminosity function with age at L^^ > 5 x 10^** erg s"' for the XB population of early-type galaxies; for older 
stellar populations, there are fewer stars > IMq, which are required to form the more luminous sources. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: individual (Centaurus A, NGC 5128) — 
X-rays: galaxies — X-rays: binaries 



1. EVITRODUCTION 

Population studies of extragalactic X-ray binaries (XBs) 
located beyond the Local Group have been made possible 
thanks to the excellent sensiti vity and spatial resolu tion of the 
Chandra X-ray observatory (iWeisskopf et al.ll2000 ). For over 
a decade, astronomers have resolved the hard X-ray spectral 
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component of galaxies at sub-arcsecon d precision into pre - 
dominantly non-nuclear point sources (ISarazinet aLll2000l) . 
Much research has focused on associations with globular clus- 
ters, which contain the m ajority of X-ray poin t sources in 
some elliptical galaxies (e.g. lAngelini et al.l2001l) . It has been 
observed that metal ri ch clusters are clea rly favored as the 
hosts of X-ray sources (iKundu et alJl2002h and t hat the dens- 
est clusters preferentially host X-ray sources (iJordan et alj 
2007). Ultraluminous X-ray sources, which display isotropic 
luminosities > lO"*^ ergs"' (Fabbiano 1989), have been 
studied in great detail by both Chandra and XMM-Newton. 
Optical observations have been used to detect count erparts 
for th ese sources in star-forming galaxies (Robert s~et al.l 
l2008h . and to demonstrate association with globular clus- 
ters in a few early-type galaxies (e.g. Maccarone et al. 2007]; 
iBrassington et al.ll2012l) . 

The X-ray luminosity function (XLF) is a well-studied 
statistical characte ristic of the X-ray s ource population for 
a given galaxy. iGrimm et al.l (120021) conducted a study 
of Galactic XBs and found that the XLF of high-mass X- 
ray binaries (HMXBs) tightly follows a power law above 
lO-'^ erg s"', while that of low-mass X-ray binaries (LMXBs) 
experiences a sharp cut-off above a few lO-'^ erg s"'. Usu- 
ally there are too few counts to derive the source spectra di- 
rectly and XLFs are produced using an assumed spectrum, 
typically a power law with F ~ 1.5-1.7, which approxi- 
mates the emission from Galactic XBs. An LMXB population 
dominat es the XLF of d iscrete sources in early-type galax- 
ies (see lFabbianoll200 6l, for review), due to their older stel- 
lar population; early-type galaxies contain relatively few high 
luminosity(> 10^^ ergs"') sources compared to late-type star- 
forming galaxies. 

The tendancy of LMXBs to dominate the discrete source 



2 



Burke et al. 



XLFs in early-type galaxies can be linked to the study of XB 
populations in late-type ga laxies, where two dist inct popula- 
tions have been observed dPrestwich et al.ll2009b . One pop- 
ulation, associated primarily with the bulge, follows an XLF 
that steepens after a few x 10^^ erg s"', while the other popu- 
lation is associated with the spiral arms and has an XLF con- 
sistent with that of HMXBs. The latter population has both 
softer colours and higher luminosities than Galactic, wind- 
driven HMXBs, most of which are accretion-powered pulsars 
with high magnetic fields, and the spectra are well-described 
by relatively cool disk blackbodies (0.1-1.1 keV), suggesting 
that they too are Roche-lobe-filling accretion driven systems. 

A luminosity break at Ljc ~ 2.5-5.0 x 10-"^ erg s"' has 
been seen in the disc rete source populat ion of numerous early- 
type galaxies (e.g., ISarazin et al.ll200i1) and is potentially the 
Eddington-limited extent of the brightest neutron star (NS) 
LMXBs. The brightest sources are super-Eddington for a 1.4 
Mq NS and are probably BH XRBs or the more massive ex- 
treme of the NS LMXB population. The luminosity break was 
examined in detail over a sample of nearby elliptical galax- 
ies by Kim & Fabbiano (2010), who found that the lack of 
sources with > 5 x 10^*^ erg s~' was more pronounced for 
older (> 5 Gyr) galaxies. The observed steepening of the 
break in the XLF suggests that young early-type XB popula- 
tions are intermediate in nature between those of star-forming 
galaxies and old early-types. 

NGC 5128 (Centaurus A) is the nearest optically luminous 
early- type galaxy, at a distance of 3.7 Mpc (F errarese et al.l 
120071) . with Mb = -2 1 . 1 . A srnall la te-type galaxy is currently 
merging with Cen A ('Graham'1979'); however, the galaxies re- 
main poorly mixed (Guillen et al. 2006). The central regions 
of the galaxy are notable for the presence of vast dust lanes 
that co ntain many active star forming regions. Rejkuba et al. 
(1201 lb compared simulated color-magnitude diagrams with 
deep ACS/HST photometry and concluded that at least 70% 
of stars in NGC 5128 formed 12 ± 1 Gyr ago and a smaller 
population of more metal rich stars formed in the last 2-4 
Gyr Therefore we expect the population of XBs to be domi- 
nated by LMXBs that possess old main sequence companions 
of M < IMq with the potential for a small component from 
LMXBs with more massive companions. 

Six 100 ks Chandra observations of NGC 5128 were 
taken as part of the Cen A Very Large Project (VLP) span - 
ning the course of 2 months in 2007 (I Jordan et al.l l200'7h . 
These observations led to investigation of the source 's strik- 
ing X -r ay features s u ch as t he AGN jet ( iHardcastl e et al.l 
(l2O07h . IWorrall et al.' ^W0§), jGoodger et all (.2010.) ), radio- 
lobe shock (Croston et al. 20® and the extended gaseous 
emission (ICraftetal. 2008). The relationship between XBs 
and globular clusters was investigated by Voss et al. (2009), 
who confirmed the presence of a low-luminosity break in the 
XLF at Lx ^ 1.5-4 X lO"*^ erg s~' and a lack of GC sources 
with Lx < 3 X lO^'' erg s"'. IVoss et all (l2009h suggested that 
this dearth of faint sources may indicate that GC LMXB com- 
panions are He rich, which possess a larger critical mass ac- 
cretion rate M„;f to become transients. Sources approaching 
the corresponding luminosity to Mc-n become unstable, and 
are therefore always above or below this luminosity, leaving 
a gap in the XLF. However, it now appears more likely that 
the low-luminosity break results from a change in the disk in- 
stability criterion at low accr etion rates as the disk sp ectrum 
peaks at longer wavelengths (Ivan Haaften et aLllMl . These 
Cen A data were also used by IZhang et aLl (1201 Ih as part of 



a large sample of GC-LMXBs from many elliptical galaxies, 
and discrepancies between the GC and field XLFs were found 
across L, - 10^^- 10^''erg s"'. 

The current era of deep X-ray observations has enabled 
more detailed s tudy of i ndividual XBs be yond the Local 
Group. Brassington et al.l (2010) and Fa bbiano et al.l (120101) 
present the results of spectral fitting the brightest sources in 
the early-type galaxies NGC 3379 and NGC 4278, down to a 
limiting unabsorbed luminosity of 1 .2 x 10^*** erg s"' . These 
samples contained 8 and 7 XBs, respectively, of which 7 were 
coincident with globular clusters (GCs). These works adopted 
a diagnostic approach to spectral fitting, inferring the state 
of a source based on fitting simple spectral models, a multi- 
colored disk blackbody and power law. The true state of 
the source can be inferred from the behavior of the absorp- 
tion parameter Nh, which simulations showed to behave in a 
characteristic way depending on the true state of the source. 
This method essentially assesses the relative contribution of 
thermal and non- thermal emission t o the spectrum, a nd i s 
fully described by iBrassington et aP (1201 Oh (see section l3l2l i. 
A further study of 18 transient-type sour ces from galaxies 
NGC 3379, NGC 4278 and NGC 4697 bv B rassington et al.1 
(120121) found a host of exotic sources, including a GC ULX 
in outflow and an unusually luminous bursting source. They 
distinguished spectral states down to ^ lO^*** erg s"'. 

In Cen A we expect to detect both neutron star (NS) and 
black hole (BH) XBs, both of which can be transients. It 
is generally thought that the transient behavior occurs due to 
the disk ionization instability mechanism (DIM), which was 
first developed for explaining the dwarf nova outbursts of cat- 
aclysmic variables (Smak 1984). The essential feature of this 
model is that the viscosity of ionized gas is larger than that 
of neutral gas. Its application to X-ray binaries is compli- 
cated by irradiation of the outer accret ion disk, which is a 
consid erably more i mport ant factor (.Dubus et al.ll2001) . Re- 
cendy, ICoriat et af] (120121) tested the DIM with a population 
of Galactic sources, finding that the critical mass accretion 
rate above which sources are persistent is lower than that pre- 
dicted when irradiation is not taken into account. Black holes 
are more massive than neutron stars, but have similar radiative 
efficiencies, and a consequence is that the outer disk temper- 
atures for black hole X-ray binaries at a given orbital period 
will be smaller than those for neutron star X-ray binaries; thus 
black hole systems are far more likely to be transient than are 
neutron star systems (King et al. 1996). 

An empirical understanding of the properties of transient 
X-ray binaries has started to develop. The size of the ac- 
cretion disk should determin e the peak outburst luminos- 
ity (e.g. 'Shahb az et al.l [T998h . and this has been borne out 
as large sampl es of such transient outbursts have devel- 
oped (Portegies Zwart et al.l (|2004 : iWu et al.l (120101) ). As 
sources change in luminosity, they fo llow loops in a hardness- 
intensity diagram dMaccarone & Co ppi 2003), indicating that 
they are changing spectral shapes hysteretically. In general, 
the spectral changes occur rapidly, with X-ray binaries in out- 
burst spending most of their time in just a few spectral states. 
Remillard & McCUntock (2006) posit that BH LMXBs pos- 
sess spectra characterized by three key spectral states. The 
thermal-dominant state, where the emission appears to be 
dominated by a ^ 1 keV multicolor disk blackbody, is es- 
sentially the same as the standard geometrically thin, opti- 
cally thick accretion disk of Shakura & Sunyaev ( 1973). Dur- 
ing the ingress and egress of outburst, the source experiences 



X-ray binaries in Cen A 



3 



a hard power-law state of F ~ 1.7, the emission likely due 
to inverse Compto n scattering in an opti cally thin, geometri- 
cally thick region (* Thorne & Pricel 11975 ). Near the peak of 
outburst, some sources also exhibit a steep power law state 
r ^ 2.5 extending to MeV energies with a significant thermal 
component also present. These spectral states are also asso- 
ciated with the rapid variability of the sources (Homan et al. 
2001b and changes in the radio jet properties (iFender et al.i 
2004). 

NS LMXBs generally show similar spectral state phe- 
nomenology to BH systems (vanderKlis 1994), but have 
some differences since the NS surface provides a boundary 
layer Nonetheless, in both cases, low/hard type spectra are 
typically seen below ^ 2% of the Eddington luminosity, ex- 
cept during the hysteretic intervals near the beginning of out- 
bursts (Maccarone 2003), and steep power law states are gen- 
erally seen only at very high luminosities, near the Edding- 
ton limit (Remillard & McClintock 2006). This spectral state 
phenomenology allows, with high quality spectra, a source to 
be classified as a candidate BH XB on the basis of its hav- 
ing a cool accretion dis k at 10^^ ergs"', as discussed by 
iWhite & Mar"shalil (IT981 . iBurke et all 120121) used the Cen A 
VLP data to show that such distinctions are now possible 
for sources outside the Local Group, presenting evidence that 
CXOU J132527.6-430023 (S14, Tabte 1) is a BH LMXB. 

The relative proximity of Cen A, coupled with the superb 
quality of these data, allows us unrivaled insight into the XBs 
of an early-type galaxy. In this work we divide sources into 
thermal and non-thermally dominant states, where appropri- 
ate, down to a luminosity of 2 x lO-'^ erg s"' - reliably mea- 
suring the spectral properties of XBs at similar luminosities 
to those found in the Local Group. 

2. DATA PREPARATION 
2.1. Source Detection and Alignment 

Each of the six 100 ks observations was analysed using 
CIAO 4.3, and was reprocessed using the chandm_repw 
script. The destreak tool was used to remove the ACIS read- 
out streak, caused by the bright Cen A nucleus. Light curves 
of each event file were produced using dmextract to check 
for background flares, which were not present. To search for 
point sources, we used a 0.5-2.0 keV event file for each ob- 
servation. We used this band because the central AGN is so 
bright (6-10 count s"' in ACIS-I) that the wings of the PSF 
contain a significant number of counts, the PSF being broader 
at higher energies for Chandra. An exposure map was cre- 
ated for each file, weighted by the typical power law spec- 
trum of an LMXB, with a photon index of F = 1.7 and ab- 
sorption column at the Gal actic value of Nh = 8.4 x 10^" cm"^ 
dPickev & Lockmanlll990 ). Inside a 5' region centered on 
the Cen A nucleus, we located point sources using the CIAO 
tool wavdetect using the spectrally weighted exposure map, 
wavelet scales of 1.0 to 16.0 in steps of \/2, a threshold sig- 
nificance of 10~* and a maximum of ^ 1 false source per 
ACIS chip. All subsequent work is within 5' of the Cen A 
nu cleus, which c orresponds to ^ the half-light radius of Cen 
A(l&aliamll979). 

The six observations were aligned by applying an appropri- 
ate x-y shift to five of the aspect solution files using the CIAO 
tool reproject_events . All observations were aligned to the 
point source positions from obsID 8490, chosen because of 
the proximity of the Cen A nucleus to the ACIS-I focus. Each 
shift was calculated from the mean offset in a and 5 obtained 



performing 2" matching between point source lists from the 
two observations. To reduce the effect of false matches on 
our offset correction, we found the mean offset in a and 5 and 
then calcuated the mean offset within ±0.5" of this mean. By 
applying 5" shifts in a and 5, subsequent 2" matching found 
^ 6 false matches between source lists. 

To allow for accurate analysis of globular cluster (GC) 
sources in our subsequent work, we utilized the well known 
LMXB -GC connection to align our observations to GC po- 
sitions (iHarris et al.ll2012h . again using 2" matching (we es- 
timate an approximate GC size of ^ 2"), we calculated the 
mean x-y shift between LMXB positions in obsID 8490 and 
the globular clusters, and applied this shift to all obsID, main- 
taining the initial alignment of the X-ray data. 

We created a merged event file using the CIAO script 
merge_all and a corresponding exposure map, as before. The 
point source list produced by a subsequent run of wavdetect 
was used as a master list of source positions that was consis- 
tent with the positions of the individual runs on each obsID, 
for each observation in which the source was detected. Cir- 
cular extraction regions were produced, centered on the chip 
position of each source, whether it is detected in an obser- 
vation or not. As was the case for the work of Burke et ai] 
(120 12l) . these regions had radius rp-j, equal to 90% of the 
7 keV extraction radius at that chip position. We then ex- 
cluded regions in the jet, radio lobe and nucleus from our 
source list. Additional region files were produced, cover- 
ing the removed read-out streak of the bright central AGN. 
Background region files were then created based on annuli 
from Irpi to 4rp7, with all source regions and the read-out 
region excluded. Subsequently we tested for source confu- 
sion, which was a significant problem in observations 7798 
and 7799, where the point source population is further off- 
axis and the point-spread function is much wider Sources 
found to have a neighbor at angular half-distance d/2 inside 
rpi were given a new extraction region radius of d/2 pro- 
vided that d/2> rp2 (where rp2 is 90% PSF radius at 2keV), 
while sources with d/2 < rp2 were declared confused (Ta- 
ble [T]i. Spectral fitting using various sized extraction regions 
for isolated, off-axis sources determined that the normalisa- 
tion found from fitting was consistent between rp2 and rpj, 
but increased systematically as the extraction radius decreased 
further 

The net counts inside each extraction region were estimated 
using dmextract to the requisite properties of the source and 
background region, followed by using aprates to calculate a 
90% confidence bound for each value or, where appropriate, a 
90% confidence upper-limit. Sources were then sorted, based 
on their highest count-flux observed in any individual obser- 
vation and names assigned based on this ranking. We present 
these results in Table [T] and emphasize that these are not the 
estimated source counts, but the net source counts present in 
a given extraction region. This indicates the quality of the 
resulting spectra we extracted for a given source. 

2.2. Source Selection 

The presence of short-term variability may indicate im- 
portant spectral changes within an X-ray source, and so we 
choose to exclude sources that vary during an observation 
from our sample. To assess the variability we made us e 
of the Gregory-Loredo algorithm (i Gregorv & Loredolll992l l. 
implemented as the CIAO tool glvary^^. This produces an 

^^ \http : // cxc.harvard.edu/ ciao/ahelp/ glvary.html\ 
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optimally-binned lightcurve for each source in each observa- 
tion, and calculates an odds ratio of each lightcurve against a 
constant count-rate, with a ratio of 10 indicating the highest 
probability of variation and being consistent with a constant 
count rate. A source is defined as showing definite variability 
when this index is > 6. The maximum ratio for each source 
is reported in Table [T] Given the potential for rapidly chang- 
ing spectra within an observation, sources with a maximum 
score of 6 or higher were removed from our sample for spec- 
tral analysis. 

At this point we also flagged sources within 20" of the 
Cen A nucleus, which require more careful extraction to deal 
with the significant flux contribution present from the nucleus. 
This contribution likely exceeds the source counts in some ob- 
servations for any sources with fewer than ^ 300 counts, de- 
pending on how far off-axis the nucleus is in the observation. 
Analysis of the sources that showed definite intra-observation 
variability and the sources within 20" of the nucleus will be 
reported in a future work. 

Table [T] summarises some basic source data; the celestial 
position, net counts in each source region for each observa- 
tion and the intra-observational variability (G - L,nax)- We also 
include our eventual source classification of persistent (P) or 
transient (T) black hole candidates (BHC), neutron star can- 
didates (NSC), foreground stars (FG), active galactic nuclei 
(AGN) and also denote S48 as a highly-magnetised NSC ((3). 

The distribution of hardness ratios showed that source S50 
was far softer than the rest of the population, with S/H^ 9 
where S is the counts from 0.5-2.0 keV and H is the counts 
from 2.0-8.0 keV. The spectru m of this source was wel l 
fit with an absorbed apec model jRavmond & Smithlll977h . 
peaking at 0.9 keV, across all observations, consistent with 
expectations from a foreground star Inspection of the data 
taken with the Inamori Magellan Areal Camera and Spectro- 
graph (IMA CS) camera on t he Magellan Observatory Baade 
telescope ( Harris et al.ll2012l) . showed that S50 is coincident 
with an object consistent with being stellar. We used ishape, 
which determines an objec t's shape pa rameters by analyti- 
cally convolving a King62 King" 196?) model with the PSF 
of the image, which is then subtracted from the input image 
of the object itself. From the residual image, the pixels are as- 
signed a weighting based on their deviation from other pixels 
at the same distance from the center of the object, and then a 
reduced is calculated. The initial parameters of the model 
are adjusted and the process is repeated with the new model 
until a minimation of the is obtained and convergence is 
reached. Our best fit with ishape indicates that our object has 
a FWHM of 0.02 pc, consistent with a star, and very different 
from the typical globula r cluster which are generally 2-4 pc in 
size ({Harris et al.ll2010h . and we do not consider it further. 
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TABLE 1 
Basic Source Properties 



Net Counts (0.5 - 8.0 keV) 



Source a S ll'il 7798 7799 7800 8489 8490 G-L^ax Type Notes 
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D. Burke etal. (2012) 
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Note. — Type column is explained in section lT2] see also the extended discussion in section|4] N: Within 20" of the Cen A nucleus, R: Source is coincident 
with read-out streak, C: Source is coincident with chip-edge, FoV: Source is outside of the field-of-view, SC: Source confused, GC: Coincident with globular 
cluster, D: Coincident with dust lane 
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3. SPECTRAL ANALYSIS 

Our primary objective is using spectral fitting of sample 
models to infer the spectral state of each source. Ideally, 
this would involve fitting all observations of a given source 
together, taking advantage of more bins to reduce the high- 
likelihood region of the parameter space. A more constrained 
absorption parameter, Nh provid es a stronger argument for 
a given spectral state (see section [J!2] |. However, it is often 
the case with joint fitting that spectra from different states of 
a source will produce a best fit that, while well constrained 
and even statistically 'acceptable', may poorly represent the 
true spectra individually. To address this issue, we designed a 
comparative test between the spectra to allow us to judge the 
extent of spectral changes for each source, and thus the appro- 
priateness of joint fits. Since we allow for changes in luminos- 
ity for a given spectral state, and only sought to determine a 
difference in the spectral shape of a source, the normalization 
of each spectrum is kept free in the joint fits. 

Spectra were extracted using the CIAO script specextmct 
for each observation that had > 100 net counts. We chose to 
fit only spectra where we had either > 150 counts in one ob- 
servation, or more than two observations with > 100 counts. 

3.L Spectral Variation 

Using the spectral fitting package XSPEC ( lArnaudlll996h . 
we performed spectral fits using an absorbed power law to 
the ungrouped spectrum, then used the covariance matrix of 
each fit in a subsequent Markov chain Monte Carlo (MCMC) 
exploration of the parameter space for each spectrum of each 
source in the sample. 



7798-8489 



3.5 - 




I I 
0.5 1.0 1.5 

Nh 

Fig. 1 . — Example MCMC output binned over Nh and F for two spectra 
of S24. We show 90% confidence regions, witli tlie overlap indicated by the 
lighter region. There is an ~ 75% probability that true combinations oi Nh 
and r are shared between the two MCMCs. 

We then tried to compare each spectrum with the (up to) 
five other spectra for a given source. For the MCMC results 



of a particular spectrum, we binned the parameter values of 
Nh and F, which describe the spectrum, into a 100 x 100 
grid; the resulting bin widths were used to extend the grid, 
if necessary, to cover the other set of parameter co-ordinates. 
This approach was intended to prevent the binning being too 
coarse or fine, the extent of the MCMC-explored region of 
the parameter space increasing as the number of counts in a 
given spectrum decreased. The number of counts in each bin, 
normalized by the total number of samples, gives the approx- 
imate 2D posterior probability density for the Nh - F space. 
The probability that the true coordinate is in an overlapping 
bin is the product of the sum of the probabilities for all the 
overlapping bins from both distributions. We adopted an a 
priori threshold of < 5% for determining if two spectra were 
too dissimilar to use in a joint fit, and the resulting spectral 
groups can be seen in our preliminary spectral fitting results 
(TableO. We show typical 90% confidence interval contours 
(fewest bins to contain 90% of all samples) of two similar 
spectra from S24 in figure [T] For a handful (~ 3) of instances 
when the behaviour of a source was not clear cut, i.e. spec- 
trum A and spectrum B were different, while both were con- 
sistent (> 5%) with spectrum C we took the more consistent 
of the two groupings, guided by the 90% confidence contours 
(i.e. a noticeable difference in shape inside the parameter 
space between two groups). 

3.2. Spectral Fitting 

After identifying sets of spectra for each source, we pro- 
ceeded to fit absorbed single-component mode ls to the 0.5-8.0 
keV spectra, following the prescription of Brassington et alj 
( 12010) , who showed that the results from such fits - partic- 
ularly the behavior of absorption parameter Nh - can give a 
strong indication of the actual spectral state of a given source. 
These results are summarised as a scheme for spectral fitting 
by means of a flow-chart in Figure 16 of Brassi ngton et al.l 
yblO), and we direct the interested reader to inspect this chart 
in tandem with our results. For thermally dominant spectra, 
where a disk blackbody component accounts for > 60% of 
the source flux, show a level of absorption significantly above 
the Galactic value N^"^ (jPickev & Lockmaniri99tf ) when fit 
with an absorbed power law, while a power law dominated 
spectrum will have an absorption less than the Galactic value 
when fit with an absorbed disk blackbody. If Nh is signifi- 
cantly above the Galactic value for both fits, then the source 
was deemed to have intrinsic absorption. If Nh was above 
A^^"' for the power law fit but significantly below N'j^"' for the 
diskbb model, then this indicates a thermal dominant spec- 
trum with some power law component also present. Con- 
versely, if Nh from fitting was zero for both models, then this 
was indicative of the source spectrum being dominated by a 
steep (F > 1 .7) power law component with a cool (kTin < 0.5 
keV) disc component also present. More specific effects are 
discussed in sectionH) 

We present results of fitting phabs x powerlaw and 
phabs X diskbb in Table [3] and indica te in the last column 
the spectral state(s) suggested by the Bra ssington et al.l (120101) 
systematic scheme. We emphasize that the results in this table 
demonstrate a phenomenological test. Sources judged to be in 
a thermal dominant or power law dominant state were subse- 
quently re-fit, imposing a lower limit on Nh at N'jf', and we 
refer the reader to Table |2]for the realistic source properties. 
The tendency for so many of the sources coincident with the 
dust lane to have an inferred 'intrinsic' absorption (i.e., Nh 
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from both absorbed power law and disk blackbody models 
was significantly above N^"'), meant that we tried to deter- 
mine the 'true' absorping column along those lines-of-sight, 
for a more reaUstic and useful comparison. 

3.3. Dust lane sources 

There is clearly a strong correlation between sources that 
possess intrinsic absorption, according to our method, and 
coincidence with the dust lane, as expected if the sources 
lie behind or within the dust. To define the spectral states 
for these sources, we obtain a value of Nh independently 
from the X-ray spectral fitting, and consider the extra ab- 
sorbing column in Cen A as well as the Galactic contribu- 
tion. We conve rted a K-band optical depth map of the dust 
lane (Kainulain en et alj|2009i) to A^^ assuming Ak ~ 0.09Av 
(IWhittetl9 92) andA^n ~ A„ x 2.1 x 10^' cm'^ (iGuver & Ozell 
I2009h . The mean Nh was then calculated using the IRAF tool 
imexarn^^ for each source coincident with the dust lanes, using 
circular apertures of radius 2", centered on the source posi- 
tion. The A^^^^ and N^^ from spectral fitting were compared 
to the mean Nh inferred from the optical depth map at each 
source position, to infer the true spectral state, if possible. 
Dust lane sources that we determined to be dominated by disk 
blackbody or power law spectra are included in Table |2] and 
in Figures|2]&|3]when the parameters could be constrain ed. A 
fuU discussion of these sources can be found in Section 14.11 

3.4. Inter-Observation Variability 

We can further our knowledge of these sources by taking the 
degree of inter-observational variability into account over the 
course of our 2 month snapshot. We calculated the net photon 
flux'' for each source region using the ciao tool aprates^^, 
after calculating the average effective exposure in the region 
using a 0.5-8.0 keV exposure map. To assess the variability, 
we calculate the fractional variability of the source, which we 
define for a set of fluxes F/ as {F,„a,x- F„,i„) / F„,tu: For sources 
that are below the detection threshold in some observations, 
the 90% upper-limit was used for f,„„,. We present these resuls 
in Figure m where we show the fractional variability against 
spectral parameter {kTi„ & F) for a given spectral state. The 
majority of sources display 10-40% variability over the time 
period spanned by the VLP. 

4. DISCUSSION 

In this section we first discuss each source coincident with 
the dust lane in turn, identifying the spectral state our method 
favors where possible. We then present the unabsorbed lumi- 
nosities as a function of key parameters of thermal and power 
law dominated states, and discuss the implications and source 
classification based on these results. We discuss possible tran- 
sient NS LMXBs in outburst and also those states for which 
we are not able to make an adequate diagnosis using simple 
models. Finally, we offer an explanation for an apparent en- 
hancement in the number of accreting BHCs found beyond 
the vicinity of the merged late-type galaxy. 

4.1. Dust lane sources 

For those sources that are coincident with the dust lane and 
where spectral fitting suggests a high Nh we infer the spec- 

http://stsdas.stsci.edu/cgi-bin/gethelp.cgi?imexamine 
" We emphasise that this does not involve assuming a spectral model 
http://cxc.harvard.edU/ciao4.2/ahelp/aprates.html 



tral state of a given source based on the mean absorption cal- 
culated from a K-band optical depth map. Due to the extra 
level of consideration given to the dust lane sources, we dis- 
cuss each source in turn below, and where a particular state 
is suggested we present it in Table [2] Our spectral diagnos- 
tic, inferring the s tate from fitting simple spectral models, is 
outlined in section [32] 

4.1.1. S4 

The brightest dust lane source, S4, is ^ 30" from the Cen A 
nucleus, coincident with a southern filament of the dust. From 
the optical depth map we calculate Nh = (3.64 ±0.07) x 
10^' cm~^, which, while above the Galactic value, is still be- 
low Nh^^(= 6.5 X 10^' cm"-). It seems likely that the spec- 
trum is dominated by a thermal component; the discrepancy 
in absorption is either a consequence of additional intrinsic 
absorption from material local to the source, or the result of 
some significant non-thermal component contributing to the 
spectrum. The true kTi„ may be somewhat cooler than that 
found from just fitting phabs(di skbb) , as suggested by the 
simulations of Br assington et aH (120101) . an effect we discuss 
in more detail for S12. Fixing Nh at the value obtained from 
the optical depth map preserves the cool disk temperature in 
spectral fitting. There is a large fall in flux from this source, 
the spectral counts falling to 29 in the final observation from 
over 1100 in the first; however, it is still above our detection 
limit. The low inner-disk temperature in the thermal state, 
combined with the high inter-observation variability, strongly 
suggests th at this is a BH L MXB transient candidate, simi- 
lar to S14 CBurkeetal .11201 2'). NS LMXBs at this luminosity, 
Lv > 10^^ erg s"', typically have a much flatter spectral shape 
over the Chandra band. 

4.1.2. S7 

For S7 the absorption map gave Nh = (6.52 ± 0.08) x 
10^' cm"^, which is below A^^'' but above N^^^ . This suggests 
that the source is behind the dust lane in a thermal dominant 
state, but with some power law component present. We are 
not able to determine the flux contribution of the power law 
component. Based on the simulations of Brassington et aP 
(i2010) we estimate the absolute systematic error ±0.5 keV 
on the temperature of the disk component. 

4.1.3. S12 

In the case of S12, flie infeiTed A^h = (16.5 ± 0.2) x 
10^' cm"^. This is significantly above the values found from 
spectral fitting, and allows for the source to be located within 
the dust lane, rather than behind it. If we assume that the 
source is behind the dust lane, i.e. has a line-of-sight absorp- 
tion column of Nh = (16.5 ± 0.2) x 10^' cm"^, then we have 
fitting results that were not produced by the simulation work 
of Bra ssington et al.l ( 120101) . where A^^** went to zero and 
A^™ was less than the Galactic value, which was the result 
obtained when the simulated spectrum was a dominant, steep 
power law component with a cool disk also present. Those 
simulations were based on a fairly low value of Nh, less than 
1 X 10^' cm"^, and so the tendency for N^^^ to go to zero 
for instances where N^'-' < N^"'"^'"^ is not surprising, as A^^^^ 
will nearly always be less than N^'^ and so the size of the 
high-likelihood region of the parameter space is small. 

To investigate our ability to recover the parameters of a 
heavily absorbed cool disk plus steep power law spectrum. 
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we carried out simulations in XSPEC. We assumed a heavily 
absorbed (Nh - 1.6 x lO^^ cm'^) cool disk (kTi„ - 0.4 keV) 
and a steep power law (F ^ 2.6). Sets of 200 simulated spec- 
tra were produced, varying the diskbb fraction of the total flux 
from 0.1 to 0.9. We then fitted phabs x diskbb to these spec- 
tra to look at the effect of varying the initial spectral parame- 
ters on those recovered by fitting a simple, single component 
spectrum. We found that the fits recovered the initial proper- 
ties of the disk when the simulated spectrum was thermally 
dominated, but as the fraction of flux taken by the power law 
component increased in significance, so the recovered value 
of Nh decreased, eventually reaching A^^^^ comparable with 
that found from fitting the spectra of S12. A consequence of 
lower A^^** as the power law contribution increased was that 
the recovered inner-disk temperature rose to beyond 1.5 keV. 
This fits neatly with the phabs x diskbb fit to S 12, which has a 
high value of kTi„ (=3.96+o^7 keV), implying an unphysically 
massive stellar BH. We take this result as an indication that 
our assumption about S12 is correct; that it is obscured fully 
behind the dust lane, with a line-of-sight Nh > 10^^ cm"^, 
and has a spectrum dominated by a power law with some 
small disk contribution. Further simulations suggest that high 
inner-disk temperatures (> 3 keV) are also recovered with a 
less steep power law component (r<1.5). This is clearly a 
complicated parameter space, the full properties of which are 
beyond the scope of this work to investigate. It is worth not- 
ing that these results indicate a significant under-estimation in 
the flux of the source by just fitting simple models; the ac- 
tual Lx is probably > 6 x 10^*^ erg s"', and could conceivably 
> 10-^^ erg s"', based on our simulations. Fitting more com- 
plex (two component) models failed to give acceptable fits 
or successfully accommodate both components, with poorly 
constrained parameters. 

S12 is analogous to S102 in NGC 3379 dBrassington et all 
12010'). for which simulations indicated that the spectra were 
dominated by a steep power law with some contribution from 
a cool disk, while a successful fit was achieved using a com- 
bined diskbb + powerlaw model yielding kTin ~ 0.14 and 
F - 1 .6 and L,- = 1 . 1 x lO^'^erg s"' . 

4.1.4. SI 3 

We calculate Nh = (7.95 ±0.13) x 10^' cm"^ in the vicin- 
ity of SI 3, which is below or consistent with the values 
obtained, but these are poorly constrained and the fits are sta- 
tistically poor with xi > 1-3. Conversely, Nh is above Nf}^^, 
which is significantly larger than Nf^"'. As was the case with 
S7, such results are indicative of a thermally dominant state 
with some non-thermal emission also present, by which we 
mean that the thermal state contributes > 60% of the total flux. 
The inner temperature in the disk of this source is less than 
that typically seen for NS systems at equivalent luminosities. 
This source displays high inter-observation variability, but re- 
mains in outburst throughout the length of our observations. 
This variability coupled with the low inner-disk temperature 
suggests that this system is a BHC LMXB, though exami- 
nation of further observations will be required to show more 
substantive evidence of the transient nature of this source. 

4.1.5. S14 

S14 is a transi ent source, analysi s of which has previously 
been reported in lBurke et al.l (12012b . A^™ was larger than the 
value inferred from the dust lane, A^^** was found to be con- 
sistent with this value, which led us to conclude the source 



was in a thermally dominant state. The cool disk at high lu- 
minosity, coupled with its transient nature, led us to conclude 
that the source is a BH LMXB candidate. Our fitting results 
for the first group suggests the presence of a power law com- 
ponent in the more luminous state. 

4.1.6. S24 

While S24 is coincident with the dust lanes, its position is 
not covered by our optical depth map. This being the case we 
cannot feasibly proceed further with our inference of spectral 
states based on the behavior of Nh, given that our first stage of 
spectral fitting is highly suggestive of increased line-of-sight 
absorption, above the Galactic value. 

4.1.7. S28 

For S28, we find Nh = (7.63 ± 0. 12) x lO^' cm"' from the 
K-band optical depth map. We fit two groups of spectra for 
this source. For the first group, consisting of spectra from ob- 
sIDs 7797 and 8490, we find a large uncertainty on Nh for 
both models, which is consistent with any location relative to 
the dust lanes at the 2a level. The second group, using spec- 
tra from obsIDs 7800 and 8489, have Nh significantly above 
j^Gaiactic^ which indicates that there is significant absorption 
along the line-of-sight to the source, which we attribute to the 
dust lane. For this group, A'^" is above the value obtained 
from the optical depth maps, while A^^^^ is consistent with 
this value, suggesting that the source is in a thermally domi- 
nant state. Assuming that the source is behind the dust lane, 
as the second set of spectral fitting results suggests, and as- 
suming a 1(7 knowledge of A'^'' = 0.51g 24 x 10^^ cm~^ and 
N^^^ = 0.44„ [5 X 10^^ cm"^, then, for the first group, the sit- 
uation is similar to that seen in SI 2, where A^™ is consistent 
with the line-of-sight value while N'^^^ is lower but non-zero. 
The higher, unrealistic (and poorly constrained), inner-disk 
temperature of 5.97 keV is consistent with this result, suggest- 
ing that the source is in a power law state during these epochs, 
with a minor contribution from a cool {kTin < 0.5 keV) disk 
component. 

4.1.8. S34 

S34 is close to the edge of the dust lane but is not coinci- 
dent with the lane itself. The high Nh in one of the spectral 
fitting groups, of which obsID 8490 provided the only spec- 
trum, prompted us to investigate the line-of-sight absorption 
using the optical depth map; we found Nh = (1.10 ±0.08) x 
10^' cm"^, consistent with the Galactic value. We are not 
able to deduce more about the state from this information; the 
source could have some intrinsic absorption not present in the 
other obsIDs. This behavior is reminiscent of the black hole 
in globular cluster RZ 2109 from NGC 4472, where the vari- 
ation in Nh is beli eved to be the re sult of a photoionizing, 
high-velocity wind (IShih et al.ll2010h . However, S34 is much 
fainter than this class of system, and the apparent increase 
in absorption could be an effect of incorrectly modelling the 
boundary layer emission, or the absorption may be a real ef- 
fect and the system is at a high inclination, undergoing a long 
period of dipping. 

4.1.9. S38 

S38 is only detected in obsID 8490. From the optical depth 
map, we infer A^// = (4.54 ±0.07) x 10^' cm"^ for this source. 
This is consistent with A^^^^ and significantly less than A^™, 
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suggesting a source in the thermally dominant state with an 
inner disk temperature ^ 0.65 keV, similar to S14, albeit at a 
lower luminosity of ^ 5.5 x 10^^ erg s"\ which points to it 
being a transient BH LMXB candidate. 

4.1.10. S46 

The X-rays from S46 would only experience an absorption 
column of A^H = (2.81 ±0.08) x 10^' cm if the source is be- 
hind the dust lane. N^"^ is less than this line-of-sight value, 
and consistent with ^Go/flcfie ^jjjjg j\jDbb j-gjj^jg j-q zero, sug- 
gesting a source in front of the dust lane, in a power law dom- 
inated hard state. 

4.2. BH & NS LMXBs 

None of the sources in our sample is confirmed to have 
L t > 4 X lO^*^ erg s"', compared to 3 sou rces apiece in 
NGC 3379 and NGC 4278 in the sample of lFabbiano et al.1 
(120101) . We caution that our sample is drawn from within 
the half-light radius, as opposed to D25, on account of the 
relative proximity of Cen A. In a ddition, Cen A is a s hghtly 
smaller galaxy with Mb = -21.1 (iFerrarese et al.ll2007h com- 
pared to -22.28 and - 22.02 for NGC 3379 and NGC 4278 
(ICappellari et al.l 120061) . and the LMXB population scales 
with galaxy mass to first order While we also exclude the 
sources that vary during the course of an observation or are 
close to the Cen A nucleus for a later study, it is unlikely 
that these sources have L^^ > 5 x lO-'^ erg s"' (Table ^ un- 
less the intra-observation variabiUty is extreme. Assuming 
the XLF foll ows an unbroken power law, with dN jdL ^ KU^ 
and a = -2 (Ki m & Fabbiand l2004). and assuming there are 
^ 40 sources between ~ 2 x lO^*^ erg s"' and ^ lO-'*' erg s"' 
(Table [D, then K ^ and there would be - 8 sources be- 
tween 10^^ erg s"' and 5 x 10^^ erg s"' and we would expect a 
single source in 5 x lO-''^ erg s"' to lO-'^ erg s"'. These values 
are consistent with our results. 

We present the thermally dominant and power law domi- 
nant spectral properties in Figures |2]&|3] respectively. In Fig- 
ure |2] we plot the unabsorbed disk luminosity against inner- 
disk temperature for the thermally dominant states (note that 
one point is representative of one ACIS spectrum, as we al- 
lowed for normalization, i.e. the flux, to be a free parame- 
ter during fitting), and we show illustrative bands of constant 
mass for lO M,^, 5M(^ and 2Mp; assuming a non-rotating com- 
pact object (iGierlinski & D one 2004), the width of the bands 
shows the variation with inclination d, from cosf? = 0.25 to 
cos^ = 0.7 5 (for which we use correction factors presented 
in IZhang et al.l 11997 ). and we assume /„,/ = 1.8. What is 
most striking is the apparent bimodal nature of the thermal 
state sources, between those that are consistent with, or to the 
right of our 2 M(r^ band, and thos e to the left of the 'SMq b. 
In Figure 6 of iBurke et al.l (120121) . we presented a population 
of Local Group thermal states, all of which lie to the left of 
our supposed 2Mq line. The fact that the bulk of our sources 
lie to the right of this line, with kTi„ > 1 keV in the range of 
^ lO-'^- lO-''^ erg s"', is very suggestive that our thermal 
state sample is dominated by NS LMXBs. 

We suggest tentatively that this result is reminiscent of the 
well-known 'mass gap' problem (Bailyn et al. 1998). The 
mass distribution of compact objects in transient systems 
strongly deviates from theoretical predictions, with a char- 
acteristic paucity between the most massive NSs (^ ^Mq) 
and the least massive BHs, while the mass distribution of 
BHs is seen to peak at ~ 8M0. If not a systematic effect 



TABLE 2 

Spectral fits using inferred dominant model (0.5-8.0 keV) 



Source N[j 


kT,„ 


Lv (0.5- 10.0 keV) 


xVdof 


Note 


10^2 


keV 


10^^ erg s"' 







DBB 



S3 


r\ nnO (19 


1.820-|4 


24.58j;75, 39.272-41 


120.7/140 




S4 


0.6500^ 




21.062;23, 37.493;;;^ 


206.9/203 


D 


S4 


0.59»;« 


0.400- 13 


6.15':™ 


10.0/14 


D 


S7 


0-5 lo6 


I.49011 


18.6|-|,27.2[-^ 


271.5/256 


D 


S8 


0.08005 




15.7li;4i,21.9l];?5 


135.7/132 




S8 






14.19|;37, 21.07212 


87.8/87 




S13 


0.400 [3 




10.57|;53,14.41]|3 


110.7/89 


D 


S13 




r\ ir^(\ 17 
^fS 
'^■■^"0.14 


'^■^'^L26'^--'-'l77 


37.5/30 


D 


S13 


0.35U;2i 




11.65f ;™ 


44.9/33 


D 


S14 


n -i-t[1 1 ■> 


nl 


7.00|;24,12.3 32;07 


72.0/71 


D 


S14 


r\ n 1 n ?S 

0.8l0;g 





18.87;^ 


14.6/27 


D 


ly 


n nQO.08 
0.09(,o, 


1.50m, 


c 7 cnO.85 
-^•^-^0.89' '•■'^0.8 


1 ly.j/ 1 jj 




S22 




1 47O.I6 
^■^'0.14 


493071 6 490.76 


63.6/67 




S26 


0.08005 


1.700-24 


6 1 30-83 7 02O.97 

0.77' 0.89 


51.6/60 


GC 


S28 


l.ioof 


2.090|7 


11.62]-76 11.917 

1.06' 1.1 


34.8/28 


D 


S31 


0.080-05 

U.Un 


1.080-15 

U. 13 


^'^"^0.53 ' '^■^0.53 


39.7/43 




S31 


0.12038 


0.980-27 

U.J 1 


3 l*^-^"^ 


13.3/12 




S32 


0.0800^ 

U.Un 


1.66S-H 


5.26!!-^„ 6.56^-0^ 

0.82' 0.98 


28.2/50 


GC 


S32 


0.080 13 

U.Uh 


1.650-55 


5.461-'' 7.672-04 

0.9/' l.oi 


16.0/16 


GC 


S34 


0.080 

0.08 


1.990-" 

0.3/ 


4.710-82 6. l9i^oi 

0.76' 0.92 


52.6/67 




S38 


600 52 


650-25 

"■"■'0.17 


5 483 68 
-^■^"1.70 


7.8/9 


D 


S44 


0.08008 

U.Un 


0. 15 


2.4l0^42,2.930^48 

0.4 ' 0.45 


29.3/42 


GC 


S45 


0.180 26 


1.650-63 

0.41 


2.860-65, 3.650^78 

O.D ' 0./ 


28.5/21 




PO 




r 








S3 




, „ I 'i 

1 75^'-^- 
'■^^11 


36.29^;f ,38.22^;^^ 


109.6/117 




S5 


0.11005 




29. 37|'^^, 36. 562-3^ 


180/214 


GC 


S8 


0.080-16 




21 11 

2.37 


23.2/34 




Sll 




^■•^^0.17 




17.3/23 


GC 


Sll 


0.0901^ 


^■■^^0.08 


12.441:37, 16.6];54 


72.3/73 


GC 


S16 


0.140-08 


^■•'^0.14 


10.05l;0^,14.26];62 


126.5/129 


GC 


S17 


1 n 


1 ^ ^0 1 7 
^■^^0.16 


IO.77J 16, I3.393-O7 


93.1/105 




S20 


0.170-!0 


1 74O.I8 
^■'^0.17 


8.430|«,11.25li 


83.5/113 




S22 




1 980 35 
^■^^^0.23 




13.4/14 




S25 


11021 

"■^^0.11 


2 110-47 

^■^^0.21 


6 24I-66 
"■■^^0.73 


14/15 




S29 


IfiO '2 


1.85q21 


S Sfi0.81 ft f;s0.87 


65.1/64 




S36 


0.200;23 


7 qq0.6I 
^■^^0.41 


4 96211 (1q2.11 


23.8/23 




S42 


13022 


1 69O 36 
^■"^0.30 


3 97O 87 5 06106 

-^•''0.79'-^-""0.99 


46.5/44 




S46 


llO''5 
"■^^0.11 


1.610-6| 


3-00°:?^4.28i:25 


28.2/27 


D 


S47 


1 l0 26 


1.860;53 


3.260:^,3.520;™ 


14.7/23 




S53 


300 32 


7 9q0.57 
^■^^0.46 


T 7irl.66 T qi 1.81 


16/19 




S55 


19059 

"■'^0.19 


1 720-83 
^■'^0.31 


3 341-56 4 362.05 
•^••^0.73' -^"0.99 


15.1/14 


GC 


S59 


llO^ 

"■'^0.11 


1 57O-67 
^■■^'0.26 


2 910-74 T ,40.91 
^•^^0.67'^^-'^0.67 


9.8/12 





Note. — Spectral fitting results with the appropriate simple model for ther- 
mally dominant (DBB) and power law dominant (PO) states, with Na forced to 
be above the Galactic value to reach a better estimate of the source luminosity. 
GC denotes that the source is coincident with globular cluster while D indicates 
that a source is in the vicinity of the dust lane. 
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(e.g. iFarr et d] 1201 It iKreidberg et all l2012i) . then this fea- 
ture in the mass distribution favours 'rapid' initial stellar col- 
lapse models of superno vae, where initial instab ihties grow on 
timescales of 10-20 ms ( Belczynski et aljr2012h . over models 
that require more prolonged instability growth (> 200 ms). 
To quantitatively test if these data are distributed bimodally in 
this plane, we first calculate a 'mass' for each source based 
on the observed peak luminosity, assuming cos 6* = 0.5. We 
then fitted^' both a single Gaussian, and then dual Gaussians 
to the observed mass distribution and from these fits calculate 
the Akaike's Information Criterion (AlC) (Akaike 1974) for 
both models. The dual Gaussian fit had means (standard de- 
viations) of 1.93 (0.76) Mq and 15.74 (8.46) Mq while the 
single Gaussian was centred on 5.6 (7.55) Mq. The AlC for 
the single Gaussian fit (=114.11) was larger than that of the 
dual Gaussian fit (=83.66), such that the latter was strongly 
favoured, with the single Gaussian 2.44 x 10"^ times as prob- 
able as the double Gaussian. These results are the first hint 
that the mass gap exists outside the Local Group. 

BH systems do not have thermal emission from a bound- 
ary layer, as one would expect from the NS LMXBs. The 
flatness of the S3 and S5 spectra, which our method sug- 
gests are power law dominated with F ^1.2, is more con- 
sistent with the spectra of LMXB NS Z-track sources than 
those of BH sources in the hard state at > 10^^ erg s"' . We 
know this based on simulating the spectra of NS LMXBs us- 
ing the model and parameters reported by Lin et al.l (120121) . 
who modeled the spectra of a bright NS LMXB, GX17H-2, us- 
ing a combination of disk blackbody, blackbody, and a power 
law component that is significant only on the Z-track hor- 
izontal branch. Regardless of whether this is a physically 
connect description or not, it is clearly a good phenomeno- 
logical description of the spectral shape. At similar count 
rates to S3 and S5, we find that the parameters of this multi- 
component model are poorly recovered by a two-component 
phabsidiskbb + bbody) fit, but that the shape of the spec- 
trum is typically well-described by a shallow power law of 
r ^ 1.0-1.4, which does not vary substantially between Z- 
track states over the Chandra bandpass. For what is judged to 
be the thermally dominant state of S3, kTi„ is higher than for 
the BH candidate systems. We therefore conclude that S3 and 
S5 are candidate Z-track sources. 

The power law states do not show a bimodality as we see 
for the thermal states, but they are consistent with Galac- 
tic NS and BH XBs, and there appears some slight favor- 
ing of steeper F at lower lu minosities . How ever, such an 
effect was also seen by Fabbiano et all (12010), whose limit- 
ing l uminosity for spec tral fitting was much higher In Figure 
9 of iBrassin^ton et al] (12010.) it is shown using 1000 count 
simulated spectra that the change in F with flux contribution 
favoured larger values from spectral fitting, particularly for 
cooler disks. It is conceivable that this effect becomes more 
pronounced when there are fewer counts. 

4.2.1. NS Transient Candidates 

Of the less luminous sources that do not appear to be BHCs 
we find that S25, S35 and S45 display extreme variability. 
S25 is only bright enough for spectral fitting in one obser- 
vation, where we determine the source to be in a power law 
dominant state. While the source is always above our detec- 
tion limit, the flux falls to ^ 20% of the initial value after 
obslD 7797, and the spectrum becomes dramatically softer 

^' Using the fanciionfitdistr, part of the MASS package in R 



by obslD 8490, with logioiS/H) ~ 1.1 and a 90% confidence 
lower-limit of ~ 0.6 in ob sID 8490. That a source displays 
quasi-super-soft behavior (iDi Stefano & Kori3 1200 3*) at low 
luminosity but F ^ 2 may indicate that S25 is a classical nova. 

S35 is only detected in obsID 8490, but in fitting simple 
models we are not able to constrain Nh such as to show that 
it is above, below or consistent with N^"' . All the parame- 
ters are poorly constrained, but taken at face value the F or 
kTi„ are not inconsistent with NS LXMBs. The variability 
demonstrated by the source would be consistent with a tran- 
sient atoll-type NS LMXB towards the peak of outburst. 

We believe that S45 may also be a transient NS LMXB in 
outburst. It displays significant temporal variability over the 
course of the six observations, and appears to be in a ther- 
mally dominant state with an inner-disc temperature kTj„ ^ 
1.3-2.3 keV, co nsistent with d isc temperatures obtained from 
Aql X-1 (e.g. Lin et al.ll2007l) . which we suggest is an ana- 
logue to this source. 

4.3. S48: A high magnetic field NS in a GC? 

S48 is a peculiar source. Coincidence with a 
spectroscopically-confirmed GC means that it is highly un- 
likely to be a background AGN or foreground star. The spec- 
trum is well-fit by a power law model of F ^ 0.7 experiencing 
negligible absorption, and a good fit cannot be achieved for 
the thermal model (Table|2]l. The luminosity of S48 makes it 
unlikely to be a cataclysmic variable, but it is similar to IGR 
J 1736 1-4441, a hard X-ray transient detected by INTEGRAL 
in the Galactic globular cluster NGC 6388 (Ferrign o et al.l 
1201 Ih , whic h is possibly a h igh-magnetic field binary such 
as GX 1h-4 (iPaul et al.l 120051) . and has a similar power law 
slope. The spectra of these sources are well described by a 
cut-off po wer law of F ^ 0.7 - 1 .0 with a high energy cut-off 
at 25 keV (iBozzo et al.ll20 1 iV The spectral shape is consis- 
tent with the compact object being a highly magnetized NS, 
the emission emanating from accretion columns that impact 
the NS surface. Young, hig hly magnetized radio pulsars have 
been observed in GCs by iBovles et al.l (1201 lb . so this is a 
plausible explanation. However, such sources ty pically show 
short-term modulation on a scale of a few hours (ICorbet et al.l 
2008); while we do not detect any variability inside of any 
observation (Table [T]!, such modulation might be too small to 
be detected. 

4.4. Uncharacterised Spectral States 

There are a handful of examples of spectra that could not 
be formally distinguished by our method (Section 13. 2t . In 
the case of S29, S32, S34 and S36, we are able to determine 
the spectral state for a subset of the observations but not for 
others that we had previously shown to have a different spec- 
tral shape. In the case of S29, obslD 7798, the source may 
be experiencing an increase in absorption or some interme- 
diate stage as it moves away from the hard state, as is the 
case with S34 obslD 8490. In obsIDs 7797, 7798 and 8490 
of S36, we find Nh consistent with the A^^"' for both mod- 
els. With the point estimate of A^^** much closer to zero, we 
could tentatively suggest that this source remains in a power 
law dominated state throughout, the spectra becoming harder 
at certain times. We note that the softer state of this source co- 
incides with its most luminous epoch, reminiscent of Galactic 
NS LMXBs such as 4U 1635-536, the spectra of which tend to 
harden as the source reaches its lowest outburst luminosities. 
In S31, there are two clear examples that seem to favor a ther- 
mal shape, but there is also clear evidence of a much harder 
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Fig. 2. — X-ray luminosity against inner-disk temperature for states deemed to be in tlie tliermal state. Diagonal bands represent constant mass for (left to right) 
10 M0, 5 Mq and 2 Mq non-rotating compact objects. The width of the bands shows an assumed inclination 6 of cos9 = 0.25 to cos9 = 0.75. 

State in ObsIDs 7798 & 8489 where it becomes more difficult 
for the fit to constrain the absorption column. For S32 obsID 
8490, Nh could not be constrained for either simple model. 



4.5. A different population? 

Spectral fitting of simple models indicates that 6 of the 8 
sources we highlight as being coincident with the dust lane 
(Table [T]i experience a greater line-of-sight absorption than 
that provided by Galactic column alone. The two sources that 
were consistent with N^"', S34 and S46, are on the edge of the 
dust lane, they are thus the worst candidates for coincidence 
and were labeled as such out of caution. Of sources not co- 



incident with the dust lane, only one source, S39, appears to 
experience excess absorption by our diagnostic. Inspection of 
the spectra of this source in more detail has led us to conclude 
that it is a background AGN, with a redshifted Fe line at ^ 4 
keV present in the spectrum. 

The good agreement between sources coincident with the 
dust lane and those that possess absorption surplus to the 
Galactic column strongly suggests that these sources experi- 
ence additional absorption from the dust lane, and validates 
our previous analysis in section |4] using a K-band optical 
depth map to estimate of the average Nh in the vicinity of 
each source (with the exception of S24, which is not within 
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Fig . 3 . — X-ray luminosity against photon index for those states determined 
to be power law dominated. 

the FoV of the map), and then infer the dominant component 
in each of the spectra. This technique assumes that the sources 
are completely behind the dust lane, and that the Nh does not 
vary substantially over a 2" radius region, while in reality the 
material is likely to be clumpy and non-uniformly distributed. 

Our four transient BH LMXB candidates are all dust lane 
sources, and are observed with kTi„ < 0.9 keV. Two further 
dust lane systems, S12 & S28, possess very cool (< 0.4 keV) 
disks with a strong power law component, mildly suggesting 
that they also possess BHs. This result indicates an enhance- 
ment of classic BHC transients inside of the merged late-type 
galaxy relative to the halo. We note that previous studies iden- 
tified the two known ULXs in the south-west quadrant of 
Cen A as tr ansient X-ray binaries that are strong BHCs (e.g., 
iGhosh et al. 2006; Sivakoff et al. 2008). Although their pres- 
ence strongly suggest that there is a population of transient 
BHCs outside of the merger region, neither source would have 
been identified as a transient based on the observations we 
consider in this paper As such, we exclude them from our 
analysis. 

Previously, hardness-intensity or color-color diagrams have 
been used to indicate which sources in a sample are heavily 
absorbed or not, but these were not ideal for giving an indica- 
tion as to whether the sources themselves were similar, as any 
further difference in hardness or color could be ascribed to the 
additional (unknown and unconstrained) absorption. These 
results from detailed spectral fitting clearly point to particular 
types of source being present in the late-type galaxy that are 
absent from the rest of Cen A. 

We propose a straightforward explanation of the lack of 
BHC transients outside of the late-type remnant. An empirical 
relation exists between orbital period and peak outburst lumi- 
nosity ( IWu et al.ll2010t) . with larger luminosities produced by 
the systems with longer orbital periods. For transients that 
contain BHs, this luminosity exceeds 10-"^ erg s"', which cor- 
responds to an orbital period of ~' 10 hours, based on the 
IWu et al, (.2010.) relation. For a main sequence companion 



filling its Roche lobe, the mass m in Mq will be related to th e 
orbital period P in hours by m '-^ O.llf jFrank et al.ll2002h . 
which means that the companion star must be be at least 1 
M0. It is likely that the donor is more evolved, particula rly in 
a BH LMXB transient (< O.ISMms, 'King et all [1991 . and 
therefore the actual main sequence lifetime of the ~ IMq 
donors is shorter The majority (70-80%) of stars in the halo 
of Cen A were formed 12 ± 1 Gyr ago, with the res t formed 
from some event 8-10 Gyr later (iReikuba et al.ll201l1) . There- 
fore the majority of stars in the halo have m < IMq, which 
accounts for the relative absence of bright BH LMXB tran- 
sients. Inside the late-type galaxy the star formation rate is 
at least 0.1 M(7) yr~', though could be as high as 1 Mq yr"' 
dMarconi et al.l uOOOl) : therefore there is a population of > 
IMq donors to feed the accreting BHs by Roche lobe over- 
flow. 

It has been proposed that the two ULXs in Cen A are BH 
LMXBs observed in the so-called ultraluminous state (e.g., 
[Gladstone et al]|2009l) at L, > lO^'^erg s"' and m quiescence 
with Lx < 10^^ erg s"' . These sources have not been observed 
in the 'classic' BH LMXB states (Remiflard & McClintock| 
120061) with L, 10^^ - 10^'' erg s"'; the non-UL outburst lumi- 
nosity of these sources is unknown and therefore they should 
not be compared directly with the dust lane BHCs until they 
have been observed below 10^' erg s"' but above the de- 
tectability threshold. In these relatively atypical sources, it 
is not yet clear how the orbital period and peak outburst lumi- 
nosity are related, which complicates direct comparison with 
the dust lane sources. 

If confirmed, the above explanation for the lack of BHC 
transients out side of the late-type re mnant also accounts for 
the results of Kim & Fab bianol ( 1201 0) for a sample of other 
nearby early-type galaxies, who found the population of 
sources with > 4 x lO^^erg s"' diminished significantly 
with early-type age, a steepening of the discrete source XLF. 
We recommend population synthesis work be carried out to 
investigate the effects of a declining population of transient 
BH LMXBs with stellar age. We also suggest a deep moni- 
toring campaign of other early-type galaxies with Chandra, to 
further investigate occurances of transient BHC XBs. 

5. CONCLUSIONS 

Our investigation into the X-ray binaries of Cen A has 
found: 

i) The population is mostly NS LMXBs, with 6 BH can- 
didates, four of which are inferred as such from the inner- 
disk temperature in the thermal state and inter-observational 
variability, and two from the inferred presence of a cool disk 
at a high luminosity. The proximity of Cen A, coupled with 
the depth of the VLP data, means that this is the only early- 
type galaxy where it is possible to perform detailed spec- 
tral fitting to sources in the classic XB luminosity range of 
lO-'^- 10^"* erg s"', similar to those sources we observe in the 
Milky Way. As a population of NS LMXBs outside of the Lo- 
cal Group, these are older analogues of sources that have been 
studied since the dawn of X-ray astronomy. We identify two 
sources that we suggest are Z track systems, based on their 
spectral shape, unabsorbed luminosity and persistent nature. 

ii) There is some evidence that the mass distribution of 
compact objects is bimodal, and there is tentative evidence of 
the so-called 'mass gap' between ^ 2.5 -SMq that has been 
observed in the population of Galactic transient LMXBs. This 
is the first time that the mass gap has been hinted at outside 
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Fig. 4.— 

the Local Group. 

iii) Besides the two known ULXs in Cen A, the only black 
hole candidate transient systems we can identify are found 
coincident with the dust lanes that arise from the merged late- 
type galaxy. These are four transient BH LMXB candidates 
observed in the thermal-dominant state. We propose that this 
is simply the result of stellar age; the older population of stars 
in the halo is unable to supply Roche lobe-filling companions 
that are required for the transient BHCs to have outburst peak 
> 10''**erg s"'. This also explains other recent results that 
show a steepening of the XLF of early-type galaxies with stel- 
lar age for L, > 4-5 x 10^'^erg s-'(Kim & Fabbiano 2010), 
as such sources would require more massive companions that 
have long since evolved off the main sequence. Similar anal- 
ysis applied to other early-type galaxies where the ages of the 



stellar populations are known, coupled with stellar population 
synthesis work, is required to investigate the robustness of this 
explanation. 
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TABLE 3 

Cen a Sources: Spectral fitting and State Identification 
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ObsIDs 
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L.I (0.5- 10.0 keV/0 




^DBB 

1022 ^^-2 




L_>- (0.5- 10.0 keV)™B 


x7dof 


State 










103' gj.g 5-1 




keV 


103' erg s"' 






S3 


7797,7800 


1 3"0* 
"•^■-'0.07 


1 25"" 


36.22;4, 38.22;^ 


109.6/117 


00" "5 


2 6l" 32 
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TABLE 3 — Continued 



Source 


ObsIDs 


NT 
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•^0.3 
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0.64 
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0.17 
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I 
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0.38 
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0.44 
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24 
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0.38 
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I 
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0.37 
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0.47 
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0.24 
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21.3/26 


I 


S42 
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n"-2 
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P1,N 


S44 
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0.30g-2^ 
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0.31 
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30.2/42 


o.oigi] 
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0.21 
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Tl 
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0.48g;g 
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Tl 
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0.11 
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0.00 
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0.36 
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PI 
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7797,8489,8490 
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^- -0.82 
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00" ™ 
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NA/NA 


M2 
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30"-32 
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00" "" 
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1 07"-2l 
^•"'0.17 


1 Q,0.37 7 70O.45 
i.5'Jq35,z.Z0(,43 
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S55 


7797,7798,7799 




1 72"**3 


S 341-55 4 ,£2.05 


15.1/14 


00" "" 

-0.82 


1 47O.53 
^•^^0.32 


7 77O.65 7 94O.9 
^•^'0.52'^-^^0.68 


17/14 


P1,N 


S59 


8489,8490 


1 1"'*2 
"•^^0.11 


1 57"-''' 


^•^^0.68'-'-^^0.81 


9.8/12 


00" "" 

"■""-0.82 


1 55"-''' 


7 070-67 7 ,0.83 
^■"'0.52'^-'0.62 


10.6/12 


P1,N 



Note. — Results from fitting absorbed power law (PO) and disk blackbody (DBB) models with all pa r amete rs free to vaiy and 95% confidence intervals. The state column denotes the spectral states that a a 
given set of spectra are consistent with, based on the classification scheme proposed bv'Brassingto n et aP ( 12010 1). T indicates a thermally dominant state, P indicates power law dominant state, M suggests that p 
there is significant, unmodelled emission from another component (> 30-40%), N suggests aresult was not-produced by the Brassin.gton et al. 12010) simulations and 1 indicates there is additional line-of-sight ^ 
absorption, when both simple models have Nh significantly above the Galactic value. The extent to which a component dominates is indicated by the associated number, i.e. 1 imphes the state is almost 
completely dominated by that component, while 2 indicates that there is some unmodelled component also present. Instances of Nj^^^ = 0.00"j|'g2 for some fits ai'e the result of not being able to consfi'ain an 
upper-limit in Xspec. We believe the source to be in the power law state at these times, with consistent with the Galactic value. In the state column, 'NA' denotes where no inference can be made regarding 
the state (see section l4~4l and for S48 shows where no fit could be achieved (i.e. x/dof > 2), which is discussed in section l431 
^ Minimal and maximal values of L, are shown for absorbed power law (PO) and disk blackbody (DBB) models 



